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Abstract: We present general results for the contact process by a method which applies to all tran- 
sitive graphs of bounded degree, including graphs of exponential growth. The model's infection 
rates are varied through a control parameter, for which two natural transition points are defined as: 
i. At, the value up to which the infection dies out exponentially fast if introduced at a single site, 
and ii. A^f, the threshold for the existence of an invariant measure with a non-vanishing density 
of infected sites. It is shown here that for all transitive graphs the two thresholds coincide. The 
method, which proceeds through partial differential inequalities for the infection density, yields 
also generally valid bounds on two related critical exponents. The main results discussed here 
were established by Bezuidenhout and Grimmett |9] in an extension to the continuous-time pro- 
cess of the discrete-time analysis of Aizenman and Barsky and of the partially similar results 
of Menshikov fl^ . The main novelty here is in the direct derivation of the partial differential 
inequalities by an argument which is formulated for the continuum. 
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1. Introduction and statement of the main results 

Since its introduction by Harris J3> the contact process has attracted interest as a model for the 
spread of "infection". The model undergoes a phase transition which is reached by varying the ratio 
of the infection rate to the healing rate, which in our notation is A : 1. The small A regime can be 
characterized by the finiteness of the "susceptibility", x{\), which is the total time lost to infection 
within the population if an infection is introduced at a single site. A duality argument allows to 
conclude that if xW < then the infection dies out even if initially the entire population was 
infected. As A approaches the edge of the regime {A : xW < oo} the susceptibility diverges and 
the contact process exhibits critical behavior with characteristics similar to those observed in models 
of statistical mechanics. 

Upon analysis, it turns out to be a generally valid statement that right past the point of divergence 
of xW ^ homogeneous contact process enters the phase at which there is a stationary measure 
with persistent infection. The technique presented below allows to establish this basic feature for 
the contact processes on the broad class of transitive graphs. This is the class of graphs which are 
invariant under the action of a symmetry group which acts transitively. Included in the collection 
are some graphs for which the contact process is known to exhibit more than one transition, in the 
sense explained below. In this generality, the basic properties of the model include the following: 

1. For A with xW < oo the probability that infection from a single site will persist in the 
population for time t decays exponentially in time. In models for which the rates for the 
direct transmission have suitable exponential decay, the probability that the infection would 
reach distance d away also exhibits, in that phase, exponential decay in the distance. 

2. At the edge of the above region x(^) Z' oo, i.e., the model exhibits criticality at the point 

At := sup{A|x(A) < oo} . (1) 

which is named for Temperley. 

3. For A > At the model is in the phase at which the infection persists. The threshold for 
the latter condition has been recognized by the term \h, for Hammersley. Thus the above 
statement amounts to the coincidence of the two points: At = A//. 

4. At the transition point, which can now be denoted simply Ac, the model exhibits critical 
behavior with characteristic exponents which in general are bounded by their 'mean-field' 
values. The bounds are realized in certain situations. 

In regard to the spread of infection, the contact process can be viewed as oriented percolation. That 
offers a helpful perspective, as the above characteristics are shared by transitive percolation models 
with or without orientation. It was in that context that the characteristics 1.-3. of the phase diagram 
were initially established for the discrete-time version of the models, in two different and indepen- 
dently derived methods, presented in the works of Menshrkov iii] and Aizenman and Barsky ll|]. 
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The argument of fT?l was limited to graphs of subexponential growth, such as Z^. The method of 
readily extends to transitive graphs and yields also additional information on the critical exponents, 
which follow through partial differential inequalities on which more is said below. However, both 
analyses were initially presented for only the discrete-time version of the models. The extension 
to the continuous-time contact process was accomplished in the work of Bezuidenhout and Grim- 
mett [9], through a detailed control of the (ID) continuum limit. Our main goal here is to present 
a direct extension of the method of Jl|] to the contact process in terms which are natural for the 
continuum, casting the argument in the generality of the transitive graphs of arbitrary growth rate. 



1.1. The model and its parameters 

We shall introduce the model in the context of transitive graphs. Before giving a formal description 
of the generator of the time evolution, let us set some notation. For Q = {V, E) a connected transitive 
graph with vertex set V and edge collection E, the contact process {At}tm. is a random time- 
dependent collection of subsets of V describing a set of infected sites. If the initial set of infected 
sites is given by B at time T, then the corresponding measure on {At}tev. is given by P^^'^) (•). Let 

( A\ - \ ^ if y e ^ 

where Jx^y is a translation invariant kernel with 

I J| = ^ Jx,y < OO. 

y 

The generator of the contact process At is formally given by 

Lf{A) = J2cix,A)[fiA . {x}) - fiA)] , (2) 

X 

where • denotes the symmetric difference operation, i.e., for x £ V: A • {x} = Au {x} when 
X ^ A and A • {x} = ^\{a;} when x G A. 

Often times the kernel is 1 when x and y are neighbors and otherwise; the finiteness of | J| im- 
plies that the vertex degree is finite. A simple description of the process in this case is as follows. At 
time t, the set of infected vertices is denoted hy At. A vertex heals independently with exponential 
rate 1, while uninfected vertices become infected at exponential rate A times the number of infected 
neighbors. 

Two significant quantities which reflect properties of the model are: 
1 . the infection density of the upper invariant measure 

^+(A) := lim p(^'^)(o G Aq) = lim p({">'^)(^o / 0) (3) 

T — >—oo T — >—oo 
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where the second equation is by duality and stationarity, 
2. the susceptibility 

X(A) := / E(W'0)(|yli|)dt <oo, (4) 
JO 

which, by Fubini, equals the expected value of the sum of the times lost to infection at the different 
sites. 

The contact process exhibits a number of different phases, depending on the control parameter A. 
Some of the thresholds of interest are defined as follows. 

Definition 1.1. 

At := sup{A : x(A) < cxo } 

\h ■■= sup{A : e+{X) =0} 

Xgn ■■= sup{A : P«°>'°)(o G ^t) — >0} 

t— »oo 

Their general relation is: 

At < Ah < Xgn ■ (5) 

Remarks: The two first transition points were already mentioned above. The third has appeared 
in the work of Grimmett and Newman 1 13] within the context of percolation models on products 
of regular trees and Euclidean lattices, where its analog is the threshold for the uniqueness of the 
infinite cluster. The above work motivated ref. 1E2I1 . where it was shown that Xh < Xgn for the 
contact process on regular trees of degree four or more; the proof was extended to all regular trees 
in 1. 161 and then more succinctly in |25]. 

The contact process can be viewed in terms of a graphical representation, whereby one traces the 
state of the infection over the 'space x time' graph, ^ x M. Healing events are represented by Poisson 
processes of intensity 1 on the lines of 1/ x M, and infection-transmission events are represented by 
Poisson processes of intensity XJx,y on {V x V) x M. For the latter, the setVxV appears as the 
collection of directed edges, and an event at {exy,t) represents a possible transmission from x to 
y at time t. The set of possible sources of infection for a site x at time t is the set of all points in 
y X M from which there is a path which does not backtrack in time, reaching (x, t) without passing 
through any healing event. We refer to this set as C{x, t). One may view it as the connected cluster 
of (x, t) in an oriented percolation model. For brevity we denote C = C(o, 0). For a more detailed 
description of the graphical representation picture and its relation to the self-duahty of the contact 
process we refer the reader to lilTT . 

The graphical representation highlights the strong relationship this process has with oriented per- 
colation, and yields the following interpretation of the two transition points: 



At = inf{A : E(IC7|) = cx)} 

Xh = inf{A : F([C| = 00) > 0}, 



(6) 
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where \C\ denotes the set's size. For a discrete set like At the size refers to the set's cardinality, 
whereas for a generic 5* C 1^ x M, such as C, we denote by |5| the total length of the set's vertical 
segments. 

By known arguments, the small-A phase has the following characteristics: 
Proposition 1.1. For any A < Aj- there exist some c < oo and r > such that 

e(W'0)(|^^|) < ce-*/". (7) 

Furthermore, if J2x Jo,x e"*"^'^' < oo for some e > 0, then also 

p({o},0) r\B^.^$ for some t > 0) < k e'^'' (8) 

for some k < oo and > 0, where CV represents the complement of a ball of radius r around 
the vertex a. 

The proof of the above proposition follows from a subadditivity property of the contact process 
and can be found in several places in the literature. For completeness we include a proof of the 
above proposition in Appendix JaI . 

1.2. Summary of the main results 

Among the key statements proven below is: 

Theorem 1.2. For any contact process on a transitive graph, with a translation-invariant infection- 
transmission rate Jx,y and a constant healing rate 1, 

\t = Xh. (9) 

As discussed above, for Q = the above result was established in 0]. The main method used 
there, based on fli'], readily extends to transitive graphs of subexponential growth (such graphs are 
amenable, though the converse is not true). It allows to conclude that at A < Xh the probability that 
the infection of one site will affect another decays exponentially in the distance (and also in time). 
However, if Q is non-amenable, e.g., a regular tree, exponential decay does not yet imply finiteness 
of x{X). Nevertheless, it is not difficult to extend the arguments given in Section 3.2 of fi*], which 
follows the approach of 1 1], to prove the above theorem for the full class of transitive graphs. 

In yj, certain non-linear differential inequalities were derived within a somewhat natural exten- 
sion of the model, for which one adds the possibility of spontaneous infection, at the rate h. Added 
insight is derived from the consideration of the model within the two parameter space of (A, h). The 
original model is then recovered through the limit /i ^ 0. In terms of the graphical representation 
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of the contact process, the spontaneous infection events are represented by a Poisson process on 
F X R with density h dt. 

It may be noted that while the extra parameter h has a very natural meaning for a contact process, 
in the original context of percolation it has appeared as a somewhat ad-hoc auxiliary "ghost field", 
whose introduction was motivated by an analogy with the external magnetic field of ferromagnetic 
Ising spin systems 121]. 

Keeping the terminology used in the percolation discussion, the events of spontaneous infection 
(points in space x time) will be referred to as green sites, and their collection denoted by G. The 
function ^+(A) which referred to the limiting density of infection starting from the 'all infected' 
state, finds its extension to /i > in the function: 

^(A,/i) =P(C(o,O)nG/0). (10) 

Following are some of the relevant properties of this extension of the model. 
Lemma 1.3. For any contact process on a transitive graph: 

i. At each h > there is a unique stationary state, to which the state of the system converges for all 
asymptotic initial condition ({S-t, —T) for T ^ oo) with the infection density given by the above 
function 9{X, h). 

ii. For h > 0, the function 9{X, h) is monotone in its arguments, continuous in X, and continuously 
differentiable in h. 

iii. In the limit h 0+, the function 6(X, h) yields the quantities which were introduced above for 
the h = model as follows: 

lime{X,h) = e+{X) (11) 

h\0 

Since the main idea is rather standard, we relegate the proof of the above lemma to Appendix (IbI . 
As can be seen there, the graphical representation provides the following useful expressions for 6 
and its derivative h) := d9{X, h)/dh is: 

e{X,h) = P(C(o,O)nG/0) = E(l-e-''l'^l) , (13) 
forA<AH: x{\h) = E (|C(o, 0)|; C7(o, 0) n G = ) = E e"'''^!) . (14) 



The graphical representation enables the derivation of partial differential inequalities which, via 
integration through the two-parameter space prove Thm ll.2l and provide also additional information 
about the behavior in the vicinity of the critical point, which can now be commonly denoted as 

Ac := At = Xh- 
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Theorem 1.4. For any transitive graph 

(i) for A > Ac.- 

6*+ (A) > Const. (A - Ac)^ (15) 

( ii) at X = Ac." 

0{Xc,h) > Const, h^'"^ . (16) 

The inequalities imply bounds for the associated critical exponents: 

/3>1, 6>2. (17) 

As explained above, these results are known already for both the discrete-time contact process ll]], 
and the continuous-time model |9]. The main novelty here is in the derivation for the continuous- 
time process of the partial differential inequalities which are discussed next. 

It should be noted that the critical ex pon ent bounds (fTTt are saturated for the contact process on 
regular trees (of degree three or more) 11241 b q). and also on when d is very large or just d > A 
and the kernel is sufficiently 'spread-out' la. 12311 . The discrete-time version of this statement was 
proven earlier through the combination of the results of 



1.3. The key differential inequalities 

The derivation of the above results proceeds through certain non-linear partial differential inequali- 
ties (PDI). The simplest of these is: 

|f<|J|/. (18) 

This relation, which for percolation was presented in Q], is basically known in the generality con- 
sidered here. It has been noted that (fTSb implies a critical exponent bound (7 < 1) which concerns 
the divergence rate for x as A /' A^: 

Next are partial differential inequalities which are similar to the PDI which were derived in Q] 
for the discrete-time contact process, in the context of percolation model, extending an earlier dif- 
ferential inequality of lloll. which has yielded a percolation analog of (fT5l with Ac interpreted as 

Theorem 1.5. For any contact process on a transitive graph, at Lebesgue almost every (A, h) G 
X M_|_ (due to the monotonicity of 9, the derivatives exist in this sense): 



(20) 
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and 

e < +(^2\'\J\9 + hx)^ + 9\ (21) 

In ref. |[l|], where the discrete version of the above theorem was established, it was envisioned 
that an extension to the continuum ought to be possible through a limiting argument, but the result 
may involve some more complicated coefficients. Nevertheless, as is shown here the inequalities 
are valid in a rather simple form, which is not that different from the discrete-time version. The 
proof proceeds through a finite-volume version of the statement, given in Theorem 13. 21 in Section|3l 
Theorem 1 1.5 1 is proved in Sectionj^ 

For the purpose of the derivation let us present some notions which are of general use when 
working with Poisson processes. 



2. A Poisson process differentiation formula 

The analysis is made clearer by recognizing a general expression for the derivatives of the probabil- 
ities of monotone events with respect to Poisson densities. It forms a continuum analog of 'Russo's 
formula' which applies in the discrete setting. 

Definition 2.1. Let "Kbe a measure space, and ^{p) a monotone functional on the space A4 of non- 
negative measures p{x) on X. A function K(x, p) on\x A4 is said to be the variational derivative 
of^{-) at X if for all finite positive continuous measures a on\ 

d 



-^{p + sa] 
as 



=0+ 



K{x) da{x) . {11) 



It is easy to see that when it exists, the variational derivative is unique. We denote it 

5^ 



K{x) . (23) 



6p{x) 

We shall now consider functionals of the form 

^(p) = Fpiu G F) (24) 

where F is an increasing event, i.e., one whose indicator function is a non-decreasing function of the 
configuration u; C X, and the subscript on P indicates that uj is distributed by the Poisson process 
with intensity measure p. 

Definition 2.2. Let F be an increasing event defined for the point process. A point x £\is said to 
be pivotal for F in the configuration oj if uj\{x} ^ F but uj U {x} e F. The set of pivotal points is 
denoted 

AF{uj) := {x:ujU{x} eF and lo\{x} ^ F} . (25) 
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Lemma 2.1. For any Poisson process, the probability of any increasing event F has a variational 
derivative given by 

= P,(F-;{xG AF}). (26) 

p{x) 

If the density p is non-atomic fp({x}) = 0/or all x E Xj then also 

5¥p{F) _ 



5p{x) 



{x e AF) . (27) 



Proof. Since the variational addition a is a continuous measure, a valid way to generate a random 
configuration distributed by the Poisson process at the density ps p + s a i?, to take the union of 
two configurations loq and a;^, drawn independently through a pair Poisson processes at intensities 
p and sa, correspondingly. By this construction, 

Pp,(F) -Ppo(F) = ¥{LOoiF,ujQ\Juis^F) 

= P (tJo ^ F, ujqUUs^ F, = 1) + (28) 
+ P (cjo ^ F, Wo U e F, > 2) 

where the first equality is due to monotonicity of F, and | • | denotes the cardinality of a set. The 
first of the two events in the last expression coincides with the event that /. luq ^ F, and //. is a 
one-point subset of AF(u;o)- The second term is dominated by [s Conditioning on loq, and 

using the explicit Poisson formula for ujg, one gets: 

P,,(F)-P,„(F) = sE(F^;a(AF(u;o))e-^"(^^('^°))) + 0(^2) 

= s j¥p{F^;x £ AF)da{x) + 0{s^). (29) 

The first claim now readily follows. 

If p is non-atomic then the probability that the site x seen on the right in eq. d26l is occupied, and 
thus F occurs, vanishes for each a-priori specified x G X. Hence the condition F"^ can be omitted 
from eq. (126 1 . which is thus reduced to eq. (l27l . □ 

One may note that some auxiliary conditions are required for an extension of the differentiation 
formula (l22l to apply also to the case where a is not a finite measure. E.g., for any set F which 
is measurable at infinity the pivotal set AF is a.s. empty, yet the probability of F need not be 
independent of p. 



M. Aizenman and P. Jung/On the Lower Phase Transition of the Contact Process 



11 



3. Derivation of the Partial Differential Inequalities 
3.1. A dictionary for the contact process 

We shall now translate Lemma ITTI to the situation at hand. Recall that the space x time picture of 
the contact process is described in terms of three independent Poisson processes describing the ran- 
dom healing events, at constant rate 1, the spontaneous infection events, and the random infection- 
transmissions. In discussing the partial derivatives of the corresponding probability, we allow the 
latter two processes to be inhomogeneous, i.e., of densities given by functions rather than constants: 
hx{t) and Xx,y{t) Jx,y The corresponding probability measure is denoted by Pa,/i- 

Of particular interest will be the event E = {C H G 7^ 0}, where C = C{o, 0) is the infecting 
cluster for a particular site (o, 0) . We apply in the natural way the terminology introduced in Defi- 
nition l2!2l and say that in a given configuration u) a site (x, t) € F x R is green pivotal if a change 
of the green set G at (x, t) will have an affect on whether u E or not. Likewise, we will say that 
an ordered bond e*^,, at time t, is bond pivotal for E if the presence of a transmission-event there 
will affect whether lo ^ E or not. 

Lemma 3.1. For the contact process on an arbitrary graph, 

Jy,x ^x,h{^yx bond pivotal for E) (30) 

P;^ t) is green pivotal for E) , (31) 



^x,hiC has exactly one green site) . (32) 

Proof. The first two assertions are direct consequences of Lemma lTTl 

For finite graphs, equation d32l can be understood from dSTT i. as we comment below. However a 
direct proof which is not limited by the finiteness condition can be obtained from the expression 
(fT3l for ^(A, h) which readily yields: 

= Pa,/i(C' has exactly one green site) , (33) 
where the last step is an explicit Poisson process relation. □ 

Remark: It is instructive to note that equation (I32t can be explained by d3]1) through the following 
argument. Let A = {lo : |C n G[ = 1} be the event that C has exactly one green site. Conditioned 
on A there is a uniquely defined site Y{uj), for which the event Y{uj) = {x,t) is characterized by: 



5e{X,h) 

5e{X,h) 
Shxit) 

and for models with h constant: 



M. Aizenman and P. Jung/On the Lower Phase Transition of the Contact Process 



12 



i. (x, t) € G, i.e., the site is an arrival point for the corresponding Poisson process, 

ii. in the configuration uj, (x, t) is a pivotal site for {C PI G 7^ 0}, i.e., for E. 

The above two statements refer to independent conditions: /. referring to the status of the site itself 
(or arbitrarily small intervals including it), and //. expressing a property of the configuration in the 
complement of this site. The probability of the former event (which is for any a-priori specified t) 
has density h with respect to dt. A simple approximation argument can be used to show that event 
//. is asymptotically independent of /. when the uncertainty interval is shrunk to a point. This yields 
the identity: 

E {I A 5{Y - (x, t))) = hF{{x, t) is green pivotal for E) , (34) 
which is to be interpreted in a distributional sense. 
Thus, 

P(|CnG| = l) = E(U) = ^ r E{IaS{Y -{x,t))) dt 

„. J —CO 



' ((x, t) is green pivotal for E) dt (35) 

5e{\,h) 



5h^{t) 



dt. 

h(-)=h 



where the last step is by dSTT i. Now, in case the total time duration of the space x time graph is 
finite, the last expression yields h^, and thus we obtain d32t . In this step we are applying the 
Definition 12.11 and Lemma 12.11 with a chosen to be the Lebesgue measure dt. This argument is, 
however, limited by the restriction in Lemma lZTl that the variational derivative a be finite. 



3.2. Differential inequalities at finite cutoffs 

As we just saw, certain technical issues need to be addressed in order to carry the analysis directly 
for an infinite graph, e.g., the decomposition of -^9 into the sum which appears at the last line of 
d35t is valid only if J2x I < ^^all circumvent this problem through finite time and space 

cutoffs. 

In order to apply arguments hke the one seen above, we let 6t,l the probability that infection is 
present at the origin o at time due to a spontaneous infection event which has occurred within the 
finite time interval (— T, 0] at some site within Vl := {x £ V : \x\ < L}. Equivalently, 9t,l is the 
infection probability at (o, 0) for the finite subgraph Ql with the vertex set Vl, in the state which 
results from A^t = 0- More generally, the infection probability in this state at (x, t) is denoted by 
Ot,l{x, t), and we let := max^^gv^ 6t,l{x, 0). 
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Due to the abundance of parameters, the dependence of the above quantities on (A, h) will occa- 
sionally be suppressed in the notation. As a step towards Theorem [731 we first derive the following 
finite-volume version. 

Theorem 3.2. On the finite graph, Qj^ x [— T, 0], the infection density introduced above 9t,l = 
^t,l(A, h), satisfies for A > and h > 0.' 

d d 

Q^^T,L < \AQtJ^ -qJ^^T,L (36) 

and 

eT,L < h-^^T,L + (2A2| J|^™f + hX) ^9t,l + [0^,T? , (37) 

This statement is proven in the rest of this section. We start with the first inequality, using the 
dictionary provided by Lemma lTTl 



Proof of d36b . Applying equation d30b : 

dt (38) 



A(.)=A 



d\ ' J-T 5Xyx{t) 

E "^y^^ J PA,/i(e^a; is vacant and is bond pivotal for £■) dt. 



y,x<=VL 

Spelling out the condition on the right-hand side we get: 

d 



0T,L = E Jy,xf^Px,hii^,t)^C-,E';Ciy,t)nG^(D)dt (39) 



9^ 

y,x€VL 



-T 



= J2 '^y,- f ^x,h{C{y,t)nGy^(D\{x,t) gC;E') Fx,h{{x,t) £C;E^)dt 

The conditional expectation is the average of the probability of the event that there is a connecting 
path from G to {y, t) in the complement of the cluster of sites which are reached from (o, 0), moving 
back in time, without visiting (x, t). Conditioning on the exact extent of this cluster, we see that the 
conditional probability is dominated by 6'!f^Jf{\, h). Thus, 

d 



g^eT,L{\h) < 0?^r(A,/i) Jy,-J_ FA,/.((x,t) eC;i?=)dt (40) 

y,x£VL 

= 9ipr{X,h)\J\ J2 [' ^x,h{{x,t)eC-E'=)dt (41) 

= (A,/i) |J| J-^T,L(A,/i), (42) 

which is the statement we wanted to show. □ 
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Inequality dTTl is a bit more involved. Start by breaking the event {C fl G / 0} into two cases: 

A = {C has exactly one green site} 
B = {|C7nG| > 2}. 

By Lemma im we have Pa,/i(^) = h-^9T,L- To estimate the probability of B, we split it further. 
Let us define a "gate" as a space x time point which, in a given configuration, needs to be visited 
by all paths which connect G to (o, 0). It is easy to see that: 

(a) The collection of gates is well ordered. 

(b) The last gate is either a green site or a vertex of a transmission bond, "or" taken in a non- 
exclusive sense (though the probability that both occur is zero.) We denote the former event 
as Bs and the latter as Bi,. 

Lemma 3.3. 

gmax o 

P«.(B.) < X k -er.. (43, 

Proof. As was discussed already, the last factor in (13.31) coincides with ¥\^}i{A). The probabilities of 
Bs and A will be compared here through the probability densities for the uniquely defined 'markers' 
for the two events. For A, that role is played by the unique green site in C, and the corresponding 
decomposition of its probability is given by equations d34l i and d35l . For Bs we note that conditioned 
on it there is a unique site W{ijj) G x [— T, 0] for which the event W{u}) = (x, t) has the 
following characteristics 

(WI) {x,t) G G 

(W2) (x, i) G C and there is no green site within the cluster of sites from which (o, 0) can be 

reached without visiting (x, t) 
(W3) there is a green site connected to (x, t) by a path in the complement of the above cluster. 

The condition (WI) has an infinitesimal probabiliy, of density h with respect to dt. Let $x^t denote 
the cluster described in (W2) and 9{x, t)^c be the probability of the event (W3) conditioned on that 
cluster. We have the following analog of equation (I34t : 

E{lBj{W-{x,t))) = hE[l[^^^r^G=<l>]0{x,t)^c) 

< hE (l[$,,nG=0] (1 - Oix, X (44) 

where we used the fact that < 9{x,t)^c < 9 max < 1 and hence 



9{x, t)^c ^ 9raax 
1 — 9{x, t)iS,c ~ 1 — 9max 



(45) 
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Now, as is easily seen, 

e(1[M)gC] l[C7nG=0] [1-^(^,0*^^]) = P((x,t) gC; CnG = 0) 

= P((x, t) is green pivotal for (46) 

Putting it together, we get the following analog of d35l 

r{Bs) = f" ^(^Bj{w-{x,t)))dt 

^ h Y^ P((x, t) is green pivotal for £■) x — 

^max X h ^ 6 (47) 

1 — Gmax dh 

which proves the lemma, through a comparison with d35l . □ 

Now we come to the trickiest estimate: 
Lemma 3.4. 

P^ ■.CR.VA M < r9\2|J|fl™«^('\ h^^hX 

dX 



Fx,h{Bb){X,h) < [2A2|J|0?^f (A,/i) + /iA1 |-^T,L(A,/i) (48) 



Proof. As in the last proof, we shall compare F (Bf,) with ^6t,l by expressing each of the quan- 
tities as integrals, with simple bounds relating the two integrands. The probability of Bf, would be 
decomposed similarly to that of Eg there, except that we shall also integrate over the specifics of the 
last event occurring at x before the time t. 

For a site (x, t), let Tx^t be the time of the last event at x preceding t, which can be either heal- 
ing, spontaneous infection, or an infection-transmission event into x. By properties of the Poisson 
distribution, for (x, t) specified: P (t — r^.t > u) = e^^^+'^+^l-^D", and conditioned on the value of 
Tx,t the probabilities of the three possibilities for the event, have the ratios 1 : h : X\J\. 

If ^x,t is as in the proof of the previous lemma, then let K^^^, K^l, and K^xt be the events that: 
G C and ^x,t n G = 0} and the last event at x preceding t is correspondingly: healing, 
spontaneous infection, or an infection-transmission event. Due to the independence of future from 
the past events, we have: 

P(iry : p(Ky : P(i^y =l:h: X\J\ (49) 

Now, if the event Bf, occurs there is a unique bond W for which the event W = e*^. is character- 
ized by the conditions: 
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(Wl) the bond e*^ is realized as an infection-transmission event, 

(W^2) there is no green site which connects to (o, 0) without visiting (x, t) 

(W^) the sites (x, t) and (y, t) are reached by a pair of disjoint paths from distinct green sites Qx 
and Qy, both in the complement of ^x,t- 

The existence of a site with the above characteristics is in fact equivalent to the event Bi,. Thus 
F {B}y) can be written, in a form similar to (l35l . as a sum of integrals of E ^ 5{W — e*^,) ^ . Splitting 

that further according to the characteristics of the last event at x preceding t, for which K^l is not 
an option, we get: 

/o 
P (k'^x}) ^ ((^' Tx^t) and (y, t) are disjointly connected to G \ K^^) • 

fc=2,3 

(50) 

In the statement that (x, r^^t) and (y, t) are disjointly connected to G, it is possible that one of the 

(2) 

paths has trivial length, e. g. the event K). /. 
For k = 2 the conditional probability in the last expression satisfies: 

P ^(x, Tx^t) and (y, t) are disjointly connected to G \ K^xVj — 
= P (^(y, t) is connected to G by a path avoiding (x, Tj^^^f) \K^J^ (51) 
y, t) is connected to G by a path avoiding (x, Tx^t) 

The first equality holds since under the condition K^. / the site (x, t) is itself green, and the second 
equality expresses the fact that the conditional probabiUty is not affected by the type of event which 
occurs at (x, Tx^t)- 

For k = ?, the condition that (x, r^^t) is infected can be met in two ways, since the site is at the 
end of an infection transmitting bond. By the van den Berg - Kesten inequality Q]; which applies 
to independent systems, the probability of the disjoint occurrence of two events is dominated by the 
product of their separate probabilities. Peeling off one of the factors, and then switching the value 
of A;, we obtain: 

P ^(x, Tx^t) and (y, t) are disjointly connected to G \ K^]^ < 

< 2 erp^f^ X IP ((y. t) is connected to G by a path avoiding (x, r^,t) I K^^}) (52) 

which 

After the above bounds are inserted in d50b . we use (l49l to change also the value of k in the factor 
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' [k^J^ appearing there. This results in: 

'{Bb) < [h + 2e^J^ \\J\]yc 

XJy^x j IP (^x^^t ; {(^1 is connected to G by a path avoiding (x, Tx^t) dt 



x,yeVL 

< [h + 2e^^if X\J\] A Jj,,:, r P (ej^ J) is bond pivotal for £; ) (53) 

= [h + 2ejpx^\\j\\\^eT,L. 

where the last equation is by d38b . □ 
Proof of d37t . Putting the above together, we have: 

St.l < f>..h(A) + Pi,»(B.) + Pj,k(Bi) 

< '.^ + '>T^^ + (2A=W''r/ + '.Al^^ (54) 



dh l-ejpi"^ dh V 11^,^' J Qx 

Collecting the first two terms and multiplying through by (1 — O^jf), one gets (l37l . □ 
This concludes the proof of Theorem l3.2l 

4. Analysis: from the PDI to the critical behavior 

We shall now extend the inequalities of Thm 13.21 to the infinite-volume (Theorem II. 5 1 . and then 
explain how they yield the main results stated in the introduction. 

Proof of Theorem 1731 By the monotonicity of the contact process, 

BTL{\h) <eT,2L{\h) < lim eT',L'{\h) = e{\h). (55) 

' T' ^L' — i-oo 

This relation permits us to simplify (linearize) the problem of passage to the limit, by replacing 
9Jp°ff in the inequalities and ( I37t by the limiting function 9. We get 

Q^OT,L{\h) < eiX,h)\J\—eT,L{X,h) (56) 

and 

Ot,l < h^^T,L + (2XMG + hX) —9t,l + e\ (57) 
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The finite-volume quantities are differentiable for /i > (in fact analytic in (A, •)). In order to take 
the limit, we shall interpret the inequalities in a weaker sense, as indicators of the corresponding 
relations for integrals of the quantities over dX dh against suitable test functions. General arguments 
permit to conclude that in this sense the inequalities remain valid also in the limit. 

More explicitely, through integration by parts d56b can be expressed as the relation of the fol- 
lowing Stieltjes integrals (each over M-|-) with positive, compactly supported, test functions g G 

Co(M+,M+) 



eT,L{X,h)dg 



dh < 



\j\eT,L{x,h)d[ge] 



dX , (58) 



where dg on the left is a Stieltjes integral at fixed h and d[g6] on the right is a Stieltjes integral at 
fixed A. By the bounded convergence theorem, as T, L — > oo, the integrals converge to those of the 
limit. Since the limiting function is also monotone in its arguments (A, h), the integration by parts 
can be reversed in the limit. The ultimate conclusion, allowed since the derivatives of monotone 
functions are locally absolutely integrable, is that the limiting inequality d20b holds in the sense of 
a relation holding at Lebesgue almost every (A, h). A similar argument permits to deduce dlTl) from 
d57b . thereby proving Theorem 1 1.51 □ 

The inequalities which are established in Theorem ll.5l are very close to what was proven in 

ref. d 

for the model's discrete-time version on Z'^. From this point on, the analysis of the PDI is identical, 
and it is covered by the general results of Lemma 4.1 and Lemma 5.1 of ref. 1 1], which yield the 
following statement (formulated here in the notation of ref. [1]). 

Proposition 4.1. Let M{(3, h) : i-^ M be a positive function which for h = is continuous from 
above and for h > is continuous, increasing in each of its arguments, and satisfies (in the a.e. 
sense): 

dM , , ^ dM 

w ^ ""'ah 

M < h'^+.^M'^^ + M^ (60) 

oh op 

with some < a < oo and some (p{P, h) and ip{(3, h) which are finite on compact subsets of 
X If there exists a value Pq for which 

limM{Po,h)/h = oo (61) 

h\0 

then for h\0 

M{Po,h) > ci/i^/(^+'') (62) 

and for (3 > Po 

M(/3,0) > C2\p-Po\'/\ (63) 

with some ci , C2 < oo. 
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Remark: Since at first glance it may appear surprising tiiat hard information about tiie critical 
behavior can be obtained from "soft " inequaUties like ( |5^ and ( 1601 . let us outUne here the heuristics 
behind Theorem l4.1l 

First, combining ( 15^ and ( 1601 . one gets: 

M < + 4>{P, h) h) M(i+") ^ + M\ (64) 

on on 

We shall apply this relation to study the h dependence at small h in the vicinity of /3o, which is 
analogous to our Ay. 

It may be noted that the inequality (164 1 does not add much information about the regime where 
M is linear in h since there h-^M ^ M, and thus already the first term on the right accounts for the 
left side. However, at Pq the dependence of M on h is singular, and may be given by a power law: 
M{Po, h) ^ h^/^ , with some 5 > 1{ the physicists convention for the corresponding exponent). For 
a shortcut, which is of course not made in the actual proof, let us allow such an assumption - taken 
in the literal sense that /iJ^M ^ \M. We now see that at /3o and h small, (l64l i holds not because of 
the first term on the right, but due to the presence of the second: 

(1 - \)M < <A(/3, h) i;i(3, h) m(^+'^) ^ + o(M) , (65) 

Dividing by M and integrating from /i = up, one gets (l62l . which leads to the interesting con- 
clusion that there is a gap in the allowed values of the exponent by which M may vanish: it either 
vanishes linearly in /i or at a slower power, 1/5, with 5 > 1 + a (in our case a = 1). 

Once it is known that for (3 > Po'- h) > c/i^/(^+"\ a similar treatment of (l60b yields for 
that regime 

1 dM 

(1 - < 0) Af — + o{M) . (66) 

Dividing by M, and integrating from Pq upward, one gets (l63l . In particular, one learns that 

M(/3, 0+) > for any P > Po \ 

The complete proof of Theorem 14.11 which does not rely on the power law assumption, can be 
obtained through the integration of the inequalities (l59b and (l60l through suitable regimes in the 
{P, h) plane, as is done in Lemmas 4.1 and 5.1 of yj. 

Proofs of Theorem ll.2l and Theorem 11.41 The two statements follow now by applying the princi- 
ple expressed in Theorem 14. II to the inequalities of Theorem II .51 with the correspondence: 
{9, A, h) 1-^ (M, P, h). For this purpose we note that by a simple estimate of the contact process 
on a graph with only one vertex, 9{X,h) > /i(l + /i)or equivalently h < 6 /{I — 9), and hence 
inequality (1211 1 can be brought to the form (l60b with ^ = 2\^\J\ + — 9) and a = 1. □ 
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5. Remarks 

1. The results presented here can be extended also to graphs which are only quasi-transitive in the 
following sense. The analysis can be adapted as long as it can be shown that for each bounded region 
in the (A, h) plane, there are < ci < C2 < oo such that 

^x,h{C{x, t) n G / 0) e [ci^(A, h), C29{X, h)] (67) 

uniformly in {x,t). In particular, the conclusions of Theorem 1 1 .21 hold under this 'weak inhomo- 
geneity' condition. 

2. Among the cases for which Theorem 11.21 applies to are the many graphs of exponential growth 
which are the subject of current research. These include hyperbolic tessellations, Cayley graphs of 
non-amenable groups, and exponentially growing amenable graphs such as the lamplighter group 
and the Diestel-Leader graph (see 1 18] and references therein.) An example for which the discrete- 
time version of Theorem 11.21 was recently applied is the thermodynamic limit of the small- world 
graphs, see fllll . 

3. The method and results presented here apply also to unoriented percolation models on transitive 
graphs similar to those considered here, i.e. ^ x M with one continuum dimension. Similar indepen- 
dence of the argument from the presence of orientation was noted in the previous related results on 
the contact process, 

4. A topic which our discussion did not address is whether in addition to the general Properties 
1.-4. it is also true, for contact process in the generality considered here, that the upper stationary 
infection density vanishes at the critical point. An equivalent formulation is that at A = Ac infec- 
tion from a single site will almost surely die out. Such a statement was established for Q = 'Z'^ 
in the celebrated work of Bezuidenhout and Grimmett |8], and their arguments can most likely be 
extended to all graphs of subexponential growth satisfying a certain homogeneity involving block 
structures. The only related results known to the authors for contact processes on graphs of expo- 
nential growth are those of ref. |20] - where the corresponding statement is proven for regular trees, 
and the corresponding statement for regular percolation on Cayley graphs of non-amenable groups, 
of ref. H. 

5. Finally, we note that Theorem 11.21 allows to sharpen a statement which was derived in iSl. As 
a step towards the proof that infection from a single site dies out almost surely at Xh, it is shown 
there, for contact processes on tree graphs, that 

exp(C(A)t) < e(W,o)(|^^|) < cexp(e(A)t) . (68) 

at some continuous ^(A), and c < co. It is not difficult to see that ^(A) > for A > Xh, and by 
Proposition 1 1 . 1 1 ^ f A ) < for A < At- Thus, Theorem ll. 21 ( A = At) allows to conclude that ^(•) 
actually changes sign at the transition point. 
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Appendix A: Exponential decay in the subcritical regime 

For completeness, we provide here a proof that throughout the regime A < At, which is character- 
ized by xW < 00, the probability that the infection, if is introduced at a single site, would persist 
for time t and/or spread over distance L decays exponentially in t and L. The proof uses generally 
known arguments. 

Proof of Proposition li.il By the additivity of the contact process and the transitivity of Q, 

E«°>'°)(|^i+,|) < E(^">'°)(|^t|)E({°>'°)(|^|) . (69) 
Subadditivity arguments permit to conclude that 

lim -logE«°>'°)(Mt|) = inf -logE(^°>'°)(MJ) = rj (70) 

t-*ca t t>0 t 

exists so that exp{r]t) < E^'^"^ '°\\At\) . Since A < At, it must be that ?? < 0. For < 5 < -?? we 
can find t so that E(^°^'°) (| At |) < exp((r? + S)t) for all t > t. Letting 1/t = -{r] + 5) and choosing 
c large enough completes the proof of Q. 

For®, we consider a process which is defined to be the contact process ignoring all healing 
events. In particular a\"^ C Fs°^ C F^°^ for all s < t. We have that the left-hand side of is 
bounded by 

p(W,o)(^ G Ft for some \x\ > r) + p(W,o)(^^ _^ 0)_ 

Coupling with a branching random walk starting from one particle at the origin gives the 
bound 

p(W,o)(^ G Ft) < e^'^pt{o,x) (71) 

for some constant K where pt (o, x) is the transition probability of a random walk on Q (note that this 
is a general bound which does not require A < Ac). A standard large deviations result which holds 
foj" Jo,x e+^l^l < 00, says that for all ^ > there is a c > such that '}2\x\>utPti^^ ^) — ce~^*. 
This together with (17 It gives an exponentially decaying bound on p({°}'*') (Fj,/„ n i?^ 7^ 0), whereas 
(0 implies that F^'^°^''^\Aj./u / 0) decays exponentially in r. □ 

Appendix B: Uniqueness of the invariant measure in the presence of spontaneous infection 

In this appendix we prove the basic regularity properties which were asserted for the process at 
h > and its relation with the standard h = version of the model. 
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Proof of Lemma ITTH /. In the presence of spontaneous infection, i.e., for /i > 0, the introduction of 
initial infection at time — T < 0, on a set S, will have negligible effect on the infection at time as 
T — > oo. To see this, we split the probability of infection at a site at time into two cases, (a) when 
it can be accounted for by a spontaneous infection event, and (b) when it can be present only due to 
the initial conditions: 

p(^'^)(oe^o) = P(C(o,o)) nGny X [-r,o] 7^ 0) + 

+p (C(o, 0)) n G n y X [-r, o] = 0; c{o, o)nB x {-t} ^ 0) (72) 

The second term is negligibly small in the limit T — > oo, since 

P {C{o, 0) n G n y X [-T, 0] = 0; C{o, 0) n S x {-T} / 0) = 

= E (e~'*l^(°'0)nGnVx[-T,o]|. (j{o,0) nSx {-T} / 0) < e"'*^ (73) 

which is obtained by first conditioning on the percolation structure, i.e., the bond variables and the 
healing events. Thus: 

lim p(^'-^) (o eAo) = P (C(o, 0) n G / 0) (74) 

= E(l-e-'^l^(°'°)l) = eiX,h). 

This implies assertion /. 

//. The monotonicity of ^(A, h) is a standard observation (and is valid also for the approximating 
functions), and the continuous differentiability of 9{X, h) in h, for /i > is an easy consequence 
(1741 . We turn our attention to the continuity of 6{\, h) in A. 

As explained above, the probability that events occurring earlier that T ago are of relevance is 
bounded by e^'^^. Restricting to times [— T, 0]: the probability that the cluster G fl y x [— T, 0] 
reaches a site with |x[ > uT can in turn be bounded by the estimates which are used in Appendix 
IXI These show that, under the assumption which is made on {Jx,y}, there exists < n such that the 
probability that the infection reaches the origin from a site at distance greater than uT is dominated 
by ce"^^. Hence, for L = uT: 

\9{X,h) -eT,L{\h)\ < e-^^ +ce-«^. (75) 

The continuity of 9t,l{X, h) in A is obvious (for a detailed argument see 1 17]). Since dTSl shows 
that 9t,l{', h) converges uniformly to 9{-, h) on it must be that 9 is continuous in A. 

The representation (l74l readily implies, via the monotone convergence theorem, that 



0(A,O+) := \iui9{\h) = P ([G(o, 0)| = oo) 

/i\0 



(76) 
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and for A < A^, 

^^^^ = i^O ^^^r = ^ i\C{o,0)\hcio,0)\^oo]) . ill) 

The event {|C(o, 0)| = 00} does not coincide with nT<o{o G ^q^''^^} however we claim that the 
difference is of probability zero. More expUcitly: 

< e{\{)+) -e+{\) (78) 

< P (|C(o, 0)1 = 00, 3r G (-00, 0) : C(o, 0) n F X {T} = 0) = 

since it is easy to see that if |C(o, 0) n F x [Ti , 0] | = 00, for some Ti < 0, then with probabihty 1: 
C{o, 0) n y X {T} / for all T < 0. □ 

The imbedding of the contact process within the its extended two-parameter version through the 
relations established in this Appendix plays a fundamental role in our analysis. In effect, it allows 
to relate the regimes of A < At and A > A^j, at /t = 0, by exploring the model along contours in 
the half plane {(A, h) eM.\ : h> 0}. 
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